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We consider a quantum field which is in a Hartle-Hawking state propagating in a spherically 
symmetric black hole space-time. We calculate the components of the stress tensor, renormalized 
using the Hadamard form of the Green's function, in the exterior region of this space-time. We then 
specialize these results to the case of the 'lukewarm' Riessner-Nordstrom-de Sitter black hole. 
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II. FORMAL CONSTRUCTION 



o 

o 

Q 



I 



> 

00 

o 

cn 



> 

X 
S3 



In a recent paper, [T] henceforth referred to as Paper I, 
we calculated the renormalized expectation value of the 
stress energy tensor operator (T'^) on the horizons of a 
spherically symmetric space-time. This calculation was 
performed using the Hadamard renormalization proce- 
dure (H 13] , which in our view provides the most direct 
and logical approach to the renormalization problem for 
practical calculations. 

Building on the seminal paper by Howard |4j, Ander- 
son, Hiscock and Samuel [S] developed a method to com- 
pute (T'^) in a general spherically symmetric space-time. 
Their method relied on the renormalization countert- 
erms, denoted by {T^^)dSi which were first calculated 
by Christensen [6 . This calculation in turn relied on 
the DeWitt series representation for the Green's func- 
tion, which is an asymptotic power series in inverse pow- 
ers of the mass of the field m. This series is ill-defined 
for the massless case and requires some severe modifica- 
tion in order to be applicable to this case. In paper we 
use the Hadamard renormalization procedure, which is 
well-defined for both massive and massless fields, to de- 
velop an alternate approach to that of [5] . In doing so we 
renormalise each term in the definition of the stress ten- 
sor independently making it very much easier to debug 
numerical calculations. 

This paper is organized as follows, in Sec. [H] we will 
outline our new method of constructing the stress tensor. 
The method leads to renormalized expressions which are 



readily amenable to numerical computation. In Sec. Ill 
we will outline the numerical method used to calculate 
these expressions. The results of this method for the luke- 



warm case are given in Sec. IV Finally our conclusions 
are presented in Sec. |Vj Throughout this paper we use 
the sign conventions of Misner, Thorne and Wheeler [7] 
and work in units in which SttG = h = c = ks = 
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In this section we present a new approach to the con- 
struction of (r'^)ren for a Spherically symmetric black 
hole space-time which has a Euclidean line element of 
the form 



ds2 = /(r)dr2 + 
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-dr^ + r^dr 
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dd(t>^. (2.1) 



The standard approach in the literature is to consider 
an expression for (T^,) y^nren in its entirety, then subtract 
off the renormalizing counterterms, (T'^)ds, as calcu- 
lated by Christensen [H], en masse to obtain a single ex- 
pression for (T'^)ren- This method was first developed 
by Candelas and Howard for a conformal scalar field in 
the Schwarzschild space-time [H [5] and then extended to 
case of scalar fields with arbitrary mass and coupling to 
the gravitational field, in general spherically symmetric 
space-times, by Anderson, Hiscock and Samuel [5]. 

We adopt an alternative approach. This approach be- 
gins with the definition of the stress tensor given in Pa- 
per I 



(2.2) 



where t'^'^[Wa] represents the coincidence limit of the 
differential operator r'"' 

T^"^ = (1 - 20gJ^V^V' + (2^ - V"'V„ 
- 2^V^V'' + 2eg^''V„V" 

+ ^(i?M-_ li^^P-)^ 1^2^M^^ (2.3) 



acting on the regular part of the Hadamard form for the 
Euclidean Green's function Ge{x,x'). Here 



(2.4) 



and gj^ is a bivector of parallel transport, which acts to 
parallel transport a vector at x' to a vector at x. 

The main idea of our method is that we consider the 
contribution from each of the derivative terms and from 
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Ge{x,x') to Eq. (2.2 1 individually. This involves calcu- yielding a finite expression which takes the form 



lating the mode sum expression for each of these quanti- 
ties, as well as the corresponding renormalising subtrac- 
tion terms. The latter quantity comes from consideration 
of the singular part of Hadamard form of G^;, which we 
will discuss in Section p^I B[ while the former may be ob- 
tained from the mode sum expression, given in Eq. (2.7) 
below, and will be the subject matter of Section |II A| This 
approach is, of course, completely equivalent to that of 
Anderson et al. however, we feels it affords a somewhat 
clearer and more tractable way of constructing (T'^)ren- 

A key issue in renormalisation, in particular when one 
cannot solve the radial equation in closed form, is finding 
a way of isolating the x — > a;' divergences in the compo- 
nents of {T^\,)unren which arc to cancel with those con- 
tained in the corresponding renormalisation subtraction 
terms. For calculations on the horizons of the space-time, 
radial separation turns out to be the most convenient 
choice of point separation. 

Indeed in paper I we made use of regularization via ra- 
dial separation to calculate (T'^)ren for a thermal state, 
on the horizons of a spherically symmetric space-time. In 
the region exterior to the black hole, excluding the im- 
mediate vicinity of the horizons, temporal separation is 
the favoured choice as angular separation generally leads 
to very complicated expressions due to derivatives of the 
Legendre function. We note here that temporal separa- 
tion only makes sense in this region, as the Killing vector 
d/dt become null on the horizons of the black hole and so 
temporal sepatation cannot be used to calculate horizon 
values. 

We will now briefiy outline the strategy employed in 
the construction of the components of {T^^)ren using 
temporal point splitting before proceeding to describe 
the details of the calculation. We follow the spirit of [5], 
but with the differences in approach outlined previously. 
Consider one component of (T'^)ren, 
[g'^" GE{x,x').^iy']ren Say (the square brackets de- 
note that the coincidence limit has been taken and we 
choose not to sum over v'), this is formally given by 
(with e(= t-t')^ 0) 

lim [g'^'^'GEix, x');..' - g^""' GEs^ng{x, a:');..'] , (2.5) 



\numeric 



[g^''''GE{x,x'),,,,]ren = [ff'^" Ge{x,x'),,,. 

Ge{^:^ );vu']analytic-i 

(2.6) 

where [g'^'^ GeIx, x').^^i]numeric consists of the coinci- 
dence limit of the mode sum expression minus the di- 
vergent sums over n contained in g'"' GEsing{x,x');i,^i , 
and [g^"^ G E{x,x').i,yi]anaiytic coutaius the finite remain- 
der terms from g^'^ GEsingix,x').,y^' . Note that we are 
free to take the coincidence limit in any direction we 
choose as W is smooth. 



A. Unrenormalized Expressions 

We begin with the mode sum expression for the un- 
renormalized Euclidean Green's function for a thermal 
state, with temperature T = k/2tt, in a spherically sym- 
metric space-time: 



Ge{x, x') — ^ F{n) cos(nK (e^)) 

oo 

X ^(2; + l)P,(cos7)p„i(r<)g„j(r>), (2.7) 



where Cr = t^t'jI'k =min(r, r'), r> =max(r, r'), F(0) = 
k/Stt^ and F{n) = k/Att'^, n > 0. Henceforth we choose 
to set r> = r. pni and qni are the independent solutions 
to the homogeneous version of this equation. 



11 (r^f^ 
dr \ dr 

= 



u'^kI l{l + l) 



(2., 



with pni defined to be the solution which is regular on 
the lower limit of the region under consideration, while 
Qni is regular at the upper limit of the region. C„; is fixed 
by the Wronskian condition 



Cnl 



Pni- 



dr 



Qnl- 



dp. 



dr 



1 



(2.9) 



where x'^ — x^ + eS^ and g'^" GEsingiXjX')-^^' repre- 
sents the renormalizing subtraction terms. Note that 
here we follow the conventions of [S], by performing all 
our derivatives in Lorentzian space, i.e. the derivatives in 
Eqn (2.2 ) are taken with respect to x = (t, r, 9, (f)). In the 



case of temporal point separation the divergences of the 
Green's function manifest themselves as divergent mode 
sums over n. In order to take this coincidence limit, we 
write the e — >■ divergences in g^'^ GEsingix, x')-^^' as 
divergent sums over n, and then bring these sums inside 
the mode sum over n. 

We are then in a position to take the coincidence limit 



Using this expression we will derive a formal expression 
(T'^)ren for the Hartle Hawking state in a spherically 
symmetric space-time. 

The first step in our computation is to calculate the 
required mode sum expressions constructed by acti on o f 
the appropriate covariant differentiations on Eq. (2.7 1. 



Firstly considering two derivatives taken at x, since g^'^ 
is diagonal it turns out that we only need to calculate 
derivatives of the form g'^'^GE^x, x')-f_i^ (here fj. is not 
summed over). For the case of a derivative taken at either 
point, g^'^ Ge^x, x')-f^,^> , the situation is not so straight- 
forward. To see this we consider the world function, a 
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(defined as half the square of the geodesic distance), for 
temporal splitting, which can be written as an expansion 



in e (see Section II B for details 
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(2.10) 



Inspection of this expression shows that even though the 
points are split in time, a has components in the (i, r) 
plane, due to the fact that the coefficients of the expan- 
sion depend on r; hence we will have cross components 
of g'^" in this plane. The components of the bivectors of 
parallel transport are readily calculated from the defining 
equation 



O"" gab'-a' — 0; [gab'] = gab, 



(2.11) 



and using the non-vanishing components of the connec- 
tion 



r'' = 



fir) 



ne = -rf{r); r.. = -r/(r) sin^ 



1. yi 



00 



cot 61; q,. 
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fir) 

Vir)' 
(2.12) 



These bivectors are then given as expansions \n e = t — t' 



0: 

„tt' 



1 2 ff2 

/ 8/ ^ 



^ r/2 Ml 
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= -9 



rt' _ ^ , -l^l I ,3 I ^ 1 



2 ^ I 48'' 48-''' 
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1 

384*^ 
1 

y.2 „;„2 



7-2 sin (0) 



384/ 

(/r + 4/YY') + o(6^) 

(2.13) 



Hence Gsix, x')-^^,' possesses two off diagonal com- 
ponents, namely g*"^ GEix,x')-tr' and g^* GEix,x')-rt'- 

Armed with these expressions we may now proceed to 
calculate the unrenormalized expressions for the required 
derivatives of the Green's function in the partial coinci- 
dence Umit (r -> r',9 9',(j) — 7> (j)'). We choose to 
denote a bitensor A{x,x') evaluated in this partial co- 
incidence limit by {A} and to drop the ix,x') notation 
for convenience. These expressions are straightforward 
to derive and so we will just list the final results: 



{Ge} ^^F{n) cos(nKe^)^(2Z-H l)p„,(r)g„,(r), 

(2.14) 



n=0 



1=0 



{g''GE;tt} = ^(") cos{nKer)Y,i^l + 1) 

n=0 1=0 

^ ( f I .dqniir) n^K^ ■ 
X -?rPniir) — 1 -^Pniir)qniir) , 



dr f 



(2.15) 



{9"'Gem'} = 

oo oo 

- 5**' Fir^W^^^ cosinKer) ^(2Z + l)p„,(r)g,u(r), 

n=0 1=0 

(2.16) 



{9''''GE;tr'} = 

- V' Y F{n)nKsm{nKer) ^(2Z + l)^^q„,(r). 



n=0 



1=0 



(2.17) 



ig^* ^ F{n)nKsin{nK€r) ^(2Z + l)p„i(r) 
«=o (=0 



dqniir) 



dr 



(2.18) 



{5'^'^ GE;rr'} = 

/ 00 00 

^rr ^ ^j-^-j cos(nKe^ )^(2Z + 1) 



'\n=0 



dpniir) dqniir) \ 
dr dr I ^ 
(2.19) 



{g''GE-ee} = {.9**Gb;^4 = 

F{n) cos{nner) Yi^l + 1) f ^p„/(r) '^'^"'^''^ 



n=0 1=0 
lil + 1) 
27-2 



-Pniir)qniir) , 



(2.20) 



00 



n=0 



00 



F{n)cos{nKer) > J(2Z + 1 ( ^-^-^^Pniir)q„iir) 



1=0 



(2.21) 



Here we have made use of the relation d/dt = id /dr. 
Note that we have not calculated a mode sum expression 
for {g^^G E-rr}- This is not required, since as will be seen 
in Section II C [g^^GE;rr]ren may be written in terms of 



the other renormalized derivatives of the Green's func- 
tion. 

It is well known that all of the sums over I above are 
divergent [8 ; this has to be a non-physical divergence 
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as the points are still separated. This superfluous diver- 
gence is just a reflection of the distributional nature of 
the Green's function. For our case, with temporal sepa- 
ration, it may be removed by subtracting off appropriate 
multiples of 6{t — r') which vanish when the points are 
separated. These subtraction terms are easily calculated 
for each of the above sums, and will be detailed in Sec- 
tion |III[ For the sake of simplicity, we will not explicitly 
incorporate them for the renormalization procedure at 
this stage; instead when we write down a mode sum over 
I, the subtraction term is included implicitly. 



B. Renormalization Subtraction Terms 

As was discussed earlier, the renormalisation subtrac- 
tion terms for a particular component are obtained by 
the action of the corresponding differential operator on 
the singular part of the Hadamard form of Ge , namely 

Gs^ng{x,x')^—+V\og{X<J). (2.22) 
(T 

Here A is a constant which is required to ensure that the 
argument of the logarithm is dimensionless. Following 
the convention of Christensen [5], we let A = 
where ^ ^ in for a massive field but is arbitrary when 
the field is massless [S]. Following [5], we choose to set 
/i = 1 for massless fields. We make use of a method due to 
Ottewill and Wardell [5j to obtain the desired subtraction 
terms, which we will now briefly outline. 

We begin by expanding a in terms of an expansion in 
powers of the coordinate separation of x and x'. Ax": 
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Ax^Ax/^ + Aa/3^Ax''Ax^Ax'^ 
+B^l}^s^x°' Ax^ Ax^ Ax^ + ... (2.23) 



where the coefficients Aap-^, Bap^s are given by 
. _ 1 



(2.24) 



B. 



abed 



a{ab^\P\cd) 



(2.25) 



The higher order terms are easily calculated by hand or 
by using Mathematica [10]. Here we have followed the 
standard convention of denoting symmetrization of in- 
dices by using brackets (e.g. (a/?)) and exclude indices 
from symmetrization by surrounding them by vertical 
bars (e.g. (a|/?|7))- 

Wardell and Ottewill [TT] have developed a Mathemat- 
ica notebook which allows one to expand both A^/^(a;, x') 
and V{x,x') in a coordinate expansion in powers of the 
coordinate separation of Ax", using Eq. (2.23), for arbi- 
trary point splitting to high order. This notebook thus 
allows one to obtain a series expansion of Gsingix, x') 



in powers of Ax", up to the required order to capture 
both the divergence and the finite remainder terms of 
Gsing{x,x') and its derivatives in the coincidence limit. 

The results of this procedure for temporal point split- 
ting are given in Appendix \K\ 



C. Renormalization 

We now proceed to obtain renormalized expressions for 
each of the derivatives of the Green's function required 
to construct the renormalized stress tensor. As outlined 
previously, we achieve this by writing the geometric x — )■ 
x' divergences in the subtraction terms as divergent mode 
sums over n and then taking the coincidence limit. To 
do this we make use of the following identities, which all 
follow from the first expression, found in [12^: 



cos{nKer) 



n=l 



= -iln(-.V) + 0(.^) 



in sin(nKeT-) 

ri=l 

oo 

K UK cos(nKeT-) 

n=l 
oo 

\k n^K^ sin(nKer) 

n=l 

oo 

n n^n^ cos(nKer) 



= - + 0(e); 



12 



+ 0(e2); 



0(e2 



6_ 

^ 120 



(2.26) 



It is at this stage of the calculation that the main 
advantage of our method becomes apparent. In previ- 
ous calculations [H [S], these identities were applied to 
the divergences contained in the Christensen subtraction 
terms, resulting in a collection of divergent sums over n. 
These sums then had to be distributed among the rele- 
vant mode sum expressions in a manner which rendered 
each sum convergent. The order of this distribution, how- 
ever, was far from obvious. In fact, in order to ensure the 
correct distribution, one had to first insert the WKB ap- 
proximations to the radial solutions into the mode sum 
expression, then perform these sums in the large n limit 
to see which counter-terms were needed for each sum. For 
our method, however, the order of distribution is imme- 
diate, as the correct divergent n sum for each component 
must come from that particular components renormaliz- 
ing counter term. 

We now employ this method to obtain a finite ex- 
pression for each of the individual derivatives of the 
Green's function required to construct the renormal- 
ized stress tensor. These calculations are repetitive, so 
will show the details of the calculation of one compo- 
nent, [Gsjren, and simply list the results for the others. 
We note here that for simplicity, we choose to perform 
our calculations for a space-time with a constant Ricci 
scalar R. Therfore we may define an effective field mass 
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m = -y/m^ + — l/6)i?. It is worth noting also that 
space-times with a cosmological constant that are solu- 
tions of Einstein's vacuum equations, possess a constant 
Ricci scalar. 



Using the identities in Eq. (2.261 we may write the 
subtraction terms Gssing hi the form 



n=l 



J nnj 6f 



2 ^"V2.2 



12 24/ 6r' 
(2.27) 



Hence we can write down an expression for [GE]ren for 
which the coincidence limit can now be readily taken: 



[GE]ren = lim < — j cos(nKe^) 



/ oo ^ 
1^(21 + l)pruqnl + J 



/ 2nK 



87r2 

(=0 

J'/ -f/2 



'o;9o; 



m' . /A/ 



In 

487r2/ 167r2 ^2 



/' 



967r2 1927r2/ 487r2r 



(2.28) 



We may now simply take the coincidence limit yielding 
the result 



[Gejren — [G e] numeric + [Ge] analytic-) 

(2.29) 



where 



\G 



E \ nu7neric 



\G 



OO 

ElanalyUc - ^^^^ ^ ^^^^ in ^^^2 J + gQ^2 
f" f 

1927r2/ 487r2r' 



H 

/ 2nK 



(2.30) 



Applying this procedure to the other components yields 
the results: 



[ff'* GE;rt']ren — [g^^ G E;tr'\ren — 0, 



(2.31) 



[5 G E\tt'\nuvaeric 



n=l 
A 



E\tt' \ analytic 



'-K^J2{2l + l)pniir)qniir) + 
A/ 



In 



/ 

Ftt' 



-r-tlK 

2 flK 



[<?- G 



E:rr' Inumeric 



47r2 



E 



4807r2/2 967r2/ 87r2/ \2k^ J Sn^f 



(2.32) 



E(2^ + i) 



dr dr 



[g Ge 



n=l ll=0 

4 



3/3 \6p 3/3 ' 2/2 ' '^'^ 



87r 



k2 //" 2/'2 m2\ Lrr' f\f , 



F I 

, J- r'T"' 



r'lanalyUc- ^^^^^2^2 96^2^3/ 3J^2 ' 
[g'^'GE-Mnumer^c = £2 E [^(^^ + 1) ^^P"' (O^r.^ (r) - - ^ 



rr 



(2.33) 



^^(2; -HI) /(; + i)poi(^)'zoi(0, 



[5^^ Ge-ob'] 



1=0 

A n ^2 



' \analytic 



Qee'K Lee' ^^j A/ \ Fee' 



14407r2/2 967r2 87r2 V2k^ / 87r2 



(2.34) 
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r„«r< 1 _ \ " 

[y ^ E-^ttlnwrneric — g^2 / y 
n=l 



^{21 + l)p„,(r)^— - — n/« + — ^{21 + l)poi(r)^^ - 





[9 G E;00\numeric = -JJ^ [[9 Ge + [.9" Gi5 

•,00' ynwmerici 

\JOr 1 - Qee^^ -^ee , / A/ \ Fee 



numeric) 



With 

Ltt, = {rV" (/" + eOm^) + 4rV" - f (rf" + 2f" - 60m') - 4rf {2f + r {r' f" + ^rf" - /")) 

+Af - 4 (45toV^ + 1) } , 
Fw = 288^^ {3rV'^ - 12rV/'^ (/" + 5m^) + Af (r^ (3rV" (/" - 30m^) - 4/'^ 

+r/' (9r/'" + 28/" - ISOm^)) + 270m + 6) + 24r/3 (2/' + r {r'f"" + 3rf"' - /")) - 24/^} . (2.37) 

Lrr' = {r'f" if" + 60m') + 4r'f" - 2r'f' {rf" + 2f" - 60m') + Arf [r {rf" + f) - 2f) 

+4/2-4(45mV + l)}, 

Frr' = 2^8^ {41rV" - 4rV/" (26/" + ISm^) + 4/^ (lOr^/" (/" + <dm') + r'f {5rf" - 24/" + 60m') 

-90m - 2) - Srf {26 f + r {ir'f"' + 15r/"' - 25/")) + 8/^} . (2.38) 

Q<^0' = {-r'f + f {r'f + 2rf + 6m'r' + 2) - 2f') 

Lee' = - {rV" + ^r'f - 2r'f {rf" + 2/" + 60m^) - 2rf (4/' + r {r'f" + 2rf" - 2/" + SOm^)) 

+4/2 + 180m^r'' + 120m'r' - 4} , 
Fee' = ^g^^ {llrV' - 8r'f (r (9/"' + rf"') + 16/") - 2r'ff' {Ur'f + lOrf + 30m'r' + lO) 

+4/2 (90m - 2) + 8f 

+4/2 (r (5r/" {r'f + em^r^ + 2) - 13r/'2 + /' (5rV" + 14r2/" + 120m2r2 + 20))) } . (2.39) 

^" = 2880rV2 i-^^^"-^" + + 2^™') + ^-^^ - ^'^^') + ^0^'^' + 2) + 2^^'-^" 

-2rV' (3r/"' + /" + 90m2) ) + 24r/3 {2f + r {r'f" + 3r/"' - /")) - 24/*} . (2.40) 

/' 

Qee = - Qee'', Lee = -Lee', (2-41) 
rf 

Fee = 2^8^ {-llrV" + 2r2//'2 (l7r2/" - 20r/' + 30m2r2 + lO) + Af {-5r'f {r'f + 6m'r' + 2) 
+13r2/'2 + rf {-5r^f" + 16r2/" + eOm^r^ - 20) + 90m'^r'' + 2) 

+8rf {30 f + r {r'f" - 6rf' - 14/") ) - 8/^ } . (2.42) 
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The construction of [g^^GE;rr]ren follows immediately from the other components. To see this we exploit the fact 
that [W] (the coincidence limit of the regular part of the Hadamard expansion for Ge) satisfies the inhomogeneous 
wave equation [3] 



{D-m^ ~^R)[W] 



-6vi 



where for a static spherically symmetric Ricci-constant space-time Vi is of the form 



1 r> nabcd 1 n nab 

^1 — abcd^ Y20 



1 



1^4 j»2 _ (^^j-, + 1) + /' (8r- - 4rV") + 4/^ + ISOtti^-* + 4} . 



1440r4 

Now since [g^^GE:rr]ren — [d^^W-rr] by definition, we have the result that 

{9 G E;rr\ren — [d G E',rr]numeric ~t~ [d G E:rr]analytici 



where 



[g^-^GE; 



rr\nu7nertc 



= -[9"GE;tt] 



numeric 



numeric + {m^+CR)[GE] numeric 



\rr\a71aiytic ^ \tt\ analytic analytic + {m^ +e,R)[GE\ analytic 



3wi 

4^' 



r 



(2.43) 



(2.44) 



(2.45) 



(2.46) 



The extension to non-Ricci constant space-times is 
straightforward and results in an extra l/nn contribu- 
tion to the sum over n from the renormalization counter- 
terms (proportional to R'). We choose not to pursue 
this extension here as we are interested in solutions to 
Einstein's equations. 



D. Formal Expressions 

In static, spherically symmetric space-times, states 
respecting the same symmetries have a stress tenor 
(T'^)renj which is diagonal. We may now write down 
formal expressions for the diagonal elements of {T^l,)ren 
in the following manner: 



in) 



analytic 



'- 1/ /ren \^ i/ /numeric 

where (not summing over v) 

(r 

ij) numeric 2(1 - Oa"" [Ge-uu' 

+ (2^ — 5) [5"" G E;aa''\numeric ~ '^SXg'^'^ G E;vi> 

•.a.oi\numeric 



(2.47) 



I numeric 



numeric 



numeric 



-[G 



E\numeric-i 



(2.48) 



and 



{"-^ I/} analytic 

2(| — Off"'' [G E;uv'\analytic + (2^ — ^)[5"" G E;aa'] analytic 
\y analytic ;aa] analytic 



analytic 



-[Ge] analytic 



2vi 
8^ 



(2.49) 



with 



1 



:R 
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The components of the Ricci tensor are given by 

f" f 

r>t r>r J _ •!_ 

n t — -n r — 2 J. ' 

R^g ^ R'*' = — 1- - 



(2.50) 



(2.51) 



If we were able to solve the radial equation ( 2.8 ) in closed 



form and perform the required mode sums we would now 
possess everything required to calculate each component 
of (r'^)ren- Unfortunately, in general, the radial equation 
must be solved numerically, adding more complexity to 
the calculation. This will be the subject matter of the 
Section Hn] 



E. Conservation Equations 

In order for the expressions we have just obtained for 
{T^^)ren to bc corrcct they must, by Wald's axioms [5], 
satisfy the conservation equations, V^(T''^)ren = 0. Us- 



ing the symmetries of the space-time and Eq. (2.12) it is 



straightforward to show that the only equation which is 
not identically satisfied is: 



8 



which, again using Eq. (2.12), reduces to 



L 

2/ 



r) Ten 



1 ren • 



(2.52) 



At this juncture we are in a position to check that the 
analytic components of the stress tensor satisfy the above 
equation. To do this we substitute the relevant expres- 
sions obtained for the analytic components obtained in 



Section II C| into Eq. (2.491 and then in turn substitute 
these into the conservation equation (2.52). Due to the 



algebraic complexity of the expression involved, this pro- 
cedure is most easily done in a Mathematica notebook 
and the result is that the right hand side of Eq. (2.52) 
vanishes. Therefore we may conclude that the analytic 
contribution to (T'^)reTi is a conserved quantity, i.e. 



analytic 



= 0. 



(2.53) 



III. 



NUMERICAL CALCULATION FOR THE 
LUKEWARM BLACK HOLE 



In this section we describe the details of the numerical 
calculations required in order to calculate both {(j>'^)ren 
and {T^^)ren in the exterior region, excluding the hori- 
zon, of a lukewarm black hole. As described in the pre- 
vious chapter, the radial equation cannot, in general, be 
solved in closed form for most space-times of interest. 
While it may be solved in special cases, for example the 
n — Q mode radial equation in Schwarzschild space-time, 
when the quantum field is massless, reduces to Legen- 
dre's equation [13j : in general, however, one has to resort 
to numerical integration to find the desired solutions. We 
begin by introducing the lukewarm black hole space-time, 
we then consider the numerical integration of the radial 
equation, before moving on the details of calculating the 
relevant mode sums. 



A. Lukewarm Black Holes 

Lukewarm black holes are a special class of Reissner- 
Nordstrom-de Sitter space-times with (Euclidean) line el- 
ement given by Eq. (2.1) with metric function 



f{r) = 1 - 



2M Q 



IT 



(3.1) 



where M, Q are the mass and charge of the black hole re- 
spectively, and A is the (positive) cosmological constant, 
with Q = M. For 4Af < a/3/A we have three distinct 
horizons, a black hole event horizon at r = r/j, an inner 
Cauchy horizon at r = r_, and a cosmological horizon at 



r = Tc, where 

r_ = 1^/3/1 ( -1 + Jl + AM^/Xji] 



(l + V^l-4A/v/A73j . 



(3.2a) 
(3.2b) 
(3.2c) 



The fourth root of / is negative and hence nonphysical. 
While the event horizon is formed by the gravitational 
potential of the black hole, the cosmological horizon is 
formed as a result of the expansion of the universe due 
to the cosmological constant [14]. An observer located 
between the two horizons is causally isolated from the 
region within the event horizon, as well as from the re- 
gion outside the cosmological horizon. 
If, as the evidence seems to suggest, the universe pos- 
sesses a cosmological constant [TS] , then it is more natural 
to consider a black hole configuration which is asymptot- 
ically de Sitter than one which sits in an asymptotically 
flat universe. Also given that de Sitter space is awash 
with radiation [13], it seems rather natural that a black 
hole in a de Sitter background would be most comfort- 
able in a final configuration in which its event horizon is 
at the same temperature as the surrounding bath. This 
state of affairs is realized in the lukewarm case and so 
the study of such a black hole configuration is well mo- 
tivated. In fact the lukewarm case has attracted much 
interest recently, as evidenced in p!6Hl8] . 
We shall confine our attention to a single exterior region 
r S [rh, rd which has a regular Euclidean section with 
topology 52 X S'2 [19]. 



B. Numerical Integration of the Radial Equation 



In this section we consider the numerical integration 



of the equation 



dS ( ^ ^dS 
— r / — 
dr V dr 



V?HiQ l{l 



(3.3) 



whose solutions will have a Wronskian satisfying 



7ll 



Pnl 



dqni 
dr 



dr 



1 

727 



(3.4) 



for a lukewarm black hole space-time as described in Sec. 



Ill A In this case, Eq. (3.3 1 has two regular singular 



points, at the event and cosmological horizons, which, 
henceforth, we will denote by and Tc, respectively. We 
apply the standard method described in Sec. V of [T7] 
to find the initial conditions for p„/ and qni about and 
Tc respectively. 

We then utilize the NDSolve algorithm in Mathematica 
to perform our integrations, with one important modifi- 
cation. NDSolve makes use of an interpolating function 



9 



to give a solution which is defined continuously on the 
range of integration. This interpolating function, how- 
ever, introduces large errors in the calculation and re- 
quires a computationally prohibative high degree of pre- 
cision in order to gain accurate results. To avoid this 
issue we force NDSolve to calculate the desired solution 
on a defined grid. We do this by splitting our integra- 
tion range into a grid, performing the integration between 
successive points and only storing the end point values. 
In particular we solve Eq. (3.3) on a evenly spaced grid 
of 1000 points to a precision of 25 decimal places. This 
procedure allows us to obtain the Wronskian of our solu- 
tions Cnl to be constant to at least 10~^^C„,/ at all grid 
points for < / < 100, < n < 7. 

Finally we choose to express all the dimensional quan- 
tities {M,Q,r,m) in units of L = -^/S/A (where A is 
the cosmological constant), which has the dimensions of 
length. 



each horizon, so we simply took the next grid point on 
each end of our region to be the horizon values calculated 
in Paper I. 



D. WKB Contribution 

We adopt the WKB approach of Howard [5, that is, 
we find an approximation to the non-linear equation that 
CniPniQui satisfies. Letting /3„;(r) = CniPni{r)qni{r) this 
equation takes the form, 



1 

2x 



'1/2 



(3.5) 



where 



C. Mode Sum Calculations Strategy 



As was discussed in Section III Al in order to ensure 
convergence of the mode sums over I one must subtract 
off multiples of the delta function, representing the di- 
vergent behaviour of the summands. These subtraction 
terms however, calculated in Section |II1D[ only render 
the mode sums 0{1/P). Therefore we would require an 
almost computationally intractable number of modes to 
achieve an accurate result. To circumvent this issue, we 
subtract off a large / approximation to the summand in- 
side the sum, thereby rendering the sum rapidly conver- 
gent. We then add back on the sum of the approximation. 
Now this sum of course also converges like 1//^, but if we 
have a closed form expression for the approximation we 
are able to transform this sum into rapidly convergent 
integrals which can be calculated with great accuracy. 

The approximation which is most commonly used for 
such calculations [H [12] is the WKB approximation, 
as it is in a relatively simple form. However, as is well 
known [5lll6j the WKB approximation suffers from issues 
with uniformity near the horizons, and as a result it fails 
to capture the correct behaviour as these horizons are 
approached, resulting in a lack of sufficient convergence 
in the numerics near the horizon. 

Therefore, the strategy we adopt is as follows: 

• For the mode sums required to calculate {T^l)ren 
we use the WKB approximation for all n modes. 

• We only use these numerical results up to a dis- 
tance from each horizon where the convergence in 
the numerical mode sums is still at a reasonable 
level. 

• We then match this last numerical point to the ex- 
act horizon values calculated in Paper I. 

In practice we found that the numerical calculations were 
of sufficient accuracy up to a distance of 1 grid point from 



(3.6) 



We desire that our WKB approximation be valid for large 
n as well as for large / so that we may cancel the renor- 
malization subtraction terms, which are given in terms of 
divergent sums over n. Hence, we look for a large x 8'P" 
proximation to i.e. a large / fixed n or a large n fixed 
I approximation. To keep track of orders it is convenient 
to replace x by x/^^ where e is an expansion parameter 
which we will ultimately set to unity at the end of our 
calculation. We then write 



(3.7) 



and expand Eqn (3.5) for small e. To balance both sides 



of this equation to lowest order we must have that 



1 

2x' 



(3.8) 



Substitution of this expression into the coefficient of 
yields an expression for which in turn gives an ex- 
pression for /32ni and so on. We may then express each 
l^ini coefficient in the expansion (3.7) in the form 



A. 



2i 



^ y2i+2j + l ■ 



(3.9) 



The Aij for the first 2 orders are given in Table [Tj the 
higher order coefficients have quite long expressions but 
are easily calculable. 

In order to calculate the mode sums required to calcu- 
late the stress tensor, we also require approximations to 
the products 



CnlPr 



dq-ai 
dr 



dpni dqni 
dr dr 



(3.10) 



10 



An n 


1 

2 






^1,0 


_i_(^Y'-4rV(r/" + 3/')- 




- I677) 


^1,1 


gLrV^^ (-3rV'' + 2r/(r/" 


+ 4/') 




^1,2 









TABLE I: WKB expansion coefficients for a spherically sym- 
metric space-time 



We can obtain the first of these expressions by exploiting 
the Wronskian condition on p„; and Qni 



Cnl 



Pnl- 



dqni 
dr 



qni - 



dpnl 

dr 



1 



and the definition of /3„; to give the identity 



CnlPr 



dqni 1 f df3,u 



dr 



dr 



1 



(3.11) 



(3.12) 



The second expression can be obtained by consideration 
of 



dr^ 



„^ dpni dqni , „ d'^q, 

^^nl— , H '^nlPnl 



dr dr 



dr^ 



Cniq-, 



dr^ 
(3.13) 



The radial equation allows us to replace the second 
derivatives of p„; and and together with the definition 
of Pni we obtain the following identity: 



^ dpni dqni _ 1 d'^jini ^ 
dr dr 2 dr^ 

'n'^K^ l{l + l) 



1/2 f\ df3,u 
2\r fjdr 



rn 



f 



Pnl. (3.14) 



We may now obtain approximations to the quantities 
(3.10) by inserting our expansion for to the required 



order, into the above identities. The error in the numer- 
ical integration of the radial equation increases with in- 
creasing I and, most severely, with increasing n; above 
n = 1 the error in the consistency of C„_; increases 
rapidly. If we choose to subtract the WKB approxima- 
tion to the order which renders the sums over / to con- 



verge like I ® and those over n like 



we need only 



sum up to a maximum of ^ = 100 and n = 7 to give 
an answer accurate to 10^^. To obtain this convergence 
we require at most a fourth order WKB approximation. 
(For many of the sums a third order approximation suf- 
fices, only when we have an P or factor is the fourth 
order approximation necessary.) We may also improve 
upon this accuracy by using the Levin u transform (see 
[20] for a detailed discussion of this transform), to speed 
up the convergence of the sums further. We now turn our 
attention to the calculation of the sums over n and I of 



the WKB approximation. We will use the transformation 

2F{i\-\) 



1=0 



Y.F{1)^ / F{\-\)d\-Tl 



1 + e 



27rA 



-d\ 



V I ^ F{l)cot(nl) , 
(3.15) 



to perform the sum over I. Here \ = I + 1/2, TZ, I 
and V denote the real, imaginary and principal parts 
respectively. We note here that Eq. (3.15) comprises of 



the standard Wastson Sommerfield formula together with 
an extra term which only contributes if F(l) takes on 
complex values for real 

To perform the required sums over / of the WKB ap- 
proximation we note that, in virtue of Eqn (3.9), each 
sum takes one of the following forms: 



2Xg,{r, n) 



/ . 1 
;=o {n^K^r^ ^\^r'^j)'2.' 



S. 



9' 



00 



2A(A2- i)/i,(r,n) 

r 

1=0 (n^K^r* + yr'^f)2 



(3.16) 



(3.17) 



where Sg and Sh represent the subtraction terms which 
are required to make each sum converge. We will consider 
sums of the first type; the same argument applies to those 
of the second form Eqn (3.17) above. Using Eq. (3.15) 
we can write 



00 



2Xgi{r, n) 

17^ 

1=0 {n^K^r^ + X^r^f)2+' 



2Xgi{r, n) 



So 



-Sn\ dX 



n 



{n^K^r'^ + X^r^f)2+' 

4A 9i{r,n) 



1 



-dX 



+1 



(3.18) 



In the terminology of Eqn (3.15), the function F{1) for 

this case is always real for real values of I so we may drop 

the extra term. Note that we drop the contribution of 

the subtraction term to the second integral above as it 

would be purely imaginary. 

We consider the real integral first. Introducing the 
1 

quantity a = nnr/f^ and changing variable from A to 
q = A/a we may recast the real integral in the form 



rf2-'0 \ (^2Ja2)i(l + g2)2 



2\n+i 



rf/'Sgj dq. 

(3.19) 
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Performing the integration gives 



It is now incumbent upon us to calculate the second 



r/2 



integral in Eq. (3.18), 



4A 



2+i 



dX 



(r2/a2)'(l-2i)- 



(3.24) 



(3.20) 

Now it is clear from the above expression that in order 
for the integral to be finite we require the subtraction 
term to satisfy the relation, 



Once more, changing independent variable from A to g 
allows us to write this integral in the form 



Ir 



n 



fhJo (i-g)i/2+^ 



dq 



(3.25) 



250 



0. 



r/1/2 
i>0 



(3.21) 



For sums of the type Eq. ( 3.17[ ) we obtain an equivalent 
expression as above with g replaced by h where 



Aaq 



Correspondingly for sums of the form Eq. (3.17), the 
above relation takes the form 



■0 (r, 2 2 



4ag(a2g2 + i) 

l + e2™9 (r2/a2)^(l + g)i/2+* 



(3.26) 
(3.27) 



^0 

Finally we arrive at the result 

„l-2i 



Inspection of Eq. (3.25) enables us to conclude that the 
integral possesses branch points at g = ±1, so we must 

1. (3-22) 

introduce a cut along the plane between [—1, 1]. 



i > 1. 



R 



2g^ 



(r/i/2)i+2* (2i~l)' 



(3.23) 



Using Eqns. (3.21) and (3.22) one can deduce the form 



of the large / subtraction terms required for each sum 
in (T'^)ren, which we denote by Sg, Scrr' and so forth. 
These are given in Table In] 



Sg 
Sow 



e 

I 1 



Ca 



1 



Son -J^('-/' + 4Av/7 + 2/ 



Sor 



+ (r - 16rC)/' + 4mV2 + SA^ + 8^ - 1) + 4(8^ + 1)/' 

(^r' if" - I6n'^) - 4/(r^(l - 4^/" 
+r(3 - 16C)/' + 4mV2 - 8X^ + 8^ + 1) + 4(8^ - l)f 



FIG. 1: The contour of integration for Ic- 

We now split the integral into 3 integrals along the con- 
tours Co, Cb and Cc, shown in Fig. [T] Thus our variable q 
to runs from — > 1 — e on Cq and from 1 + e ^ 00 on Ce- 
lt is clear that the integral along Cc is purely imaginary 
and so does not contribute to the final answer. Both the 
integral along Ca and along Cb will have a contribution 
which will be divergent in the e — > limit. In Appendix 
|B]we show that these contributions cancel out leaving a 
answer which is finite in the e — > limit, given by: 



E 



TABLE II: Large I subtraction terms 



(-i)^gF'(i) 

r/2 + 
/o 



(3.28) 
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This result is equally valid if g is replaced by h defined 
by Eq. KTf\. 



Now while Eq. (3.28) is a perfectly finite expression, 
in practice we must calculate the integral numerically, 
which will run into accuracy errors near q = \ despite 
being finite there. This is due to the divergent nature 
of the denominator in the integral, which for large i will 
diverge rapidly, so that inevitable round-off errors in the 
denominator will lead to potentially large errors in the 
numerical integration. To avoid these errors we integrate 
by parts i — 1 times (see Appendix |b] for details) leading 
to the following expression for Ic'- 



Ic = 



r/2 



.9.(0) 



2)(* 2' 



+2 



(-1)' 



Jo 



5f'(0) 



(3.29) 



which may now be easily calculated numerically. Again, 
the corresponding result, when g is replaced by h, is also 

valid^ Finally we note that for the n = mode, from 

Eq. (3.6 1, X = (^ + l/2)-\/r2/, therefore each sum over I 
is of the simple form: 



i/(r)^(Z + l/2) 



-2j 



(3.30) 



1=0 



where j is an integer > and H{r) is some function of 
r. Due to their simplicity, these sums may be directly 
evaluated with out recourse to the Watson-Sommerfeld 
method, using P^IH] 



+ = -(1 - 22")C(2n) 

(1 - 22")22""i 
" " (-l)"-i(2n!) 



B2n- 



(3.31) 



where C is Riemann's Zeta function and B2n are the 
Bernoulli numbers. 

At this juncture it is useful to take stock of our progress 
thus far. We are now in a position to calculate the sum 
over I and n of the full radial solutions minus their WKB 
approximations to a reasonable accuracy (~ 10~^). We 
may also perform the sum over / of the WKB approxima- 
tion using the integration method outlined above. There- 
fore all that remains for us to calculate in this section, is 
the sum over n of these WKB integrals. We use the Levin 
u transform to perform this sum over n. We calculate the 
first 20 n values and then apply the Levin u transform 
to these values to give a final value for the sum. 

We now have all tools required to calculate both 
{4>'^)ren and (T'^)ren numerically in the exterior region. 



excluding the immediate vicinity of the horizons. Plots 
of the results combined with the exact hoirzon values 
calculated in Paper 1, can be found in the next section. 



IV. RESULTS 

A. Plots of (f, 

In Fig. [2] we plot the components of the stress tensor, 
calculated on a grid of 1000 points, for a massless con- 
formally coupled scalar field in the exterior region of a 
lukewarm black hole with M = Q = O.IL. We note here 
that due to spherical symmetry (Tg^)ren = {T^)rem so 
we do not include a plot of the latter. {Tj^^)ren is plotted 

Components of 

(T/>„„ 
0.4 h„ 





0.2 



0.4 



0.6 



Cosmological 
Horizon 



FIG. 2: A plot of the non-zero elements of {T^^}ren in the 
entire exterior region. 

on the first 100 points on the grid plus the event horizon 
value in Fig. g while (f/)^e„ is shown on the last 100 
viable grid points plus the cosmological horizon value in 
Fig. [i] In Figs. |3]and|4]we see the equality of {Tj.'^)ren 



Components of 

(r„">,.e„ 



-(V)„ 
- (r/),. 




FIG. 3: A plot of the non-zero elements of {T^^)ren in the 
region of the event horizon. 

and {T^^)ren on both horizons as derived in Paper 1. 
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Components of 



Cosmological Horizon 

FIG. 4: A plot of the non-zero elements of {T^^)ren in the 
region of the cosmological horizon. 



Components of the 
Conservation Equatio 




C_2 
C 



Cosmoiogicai Horizon 



B. Conservation Equations 

In the Section [lIE| we demonstrated that the analytical 
contribution to {TDren is a conserved quantity. We are 
now in a position to perform the same analysis on the 
numerical contribution for a lukewarm black hole space- 
time. The conservation equation is given by 

2 - 

/' " "t 2 

+ '^{{'-^ r)ren ~ {T f}ren) ~ -^{T g)ren- 

(4.1) 

As we have performed all of our numerical calculations 
on a grid we see that in order to construct the quan- 
tity {T^^)ren,r, wc must employ a numerical derivative 
scheme. We choose to implement a finite difference 
scheme to calculate the derivative. Mathematica has a in- 
built finite difference algorithm which we employ. Unfor- 
tunately, calculating (T''"^)ren,r using this method is prone 
to error, in particular near the horizon where the rate of 
change is large. These errors introduce a lot of numerical 
noise into the conservation equation. However, a more 
informative measure, to that of the conservation equa- 
tion alone, is the magnitude of the conservation equation 
(C) relative to that of its constituent parts, which we 
denote as 



1 r/r^n^r: 




(4.2) 

In Fig. [5] we see a plot of the conservation equation and 
its components between the event and cosmological hori- 
zon for a massless scalar field on lukewarm black space- 
time with parameter values M = Q = O.IL. In Fig. [6] 
we zoom into a mid section of the exterior region, again 
for a massless scalar field on lukewarm black space-time 
with the same black hole parameter values. These plots 
show that the conservation equation is small relative to 



FIG. 5: A plot of the components of the conservation equa- 
tion over the entire exterior region. 




FIG. 6: A plot of the components of the conservation equa- 
tion over a mid-section of the exterior region. 



its constituent components, Ci and C2, almost every- 
where with the exception of the region of the point where 
Ci — C2 — and the region of the cosmological horizon. 

In the neighborhood of the cosmological horizon, the 
magnitude of the quantitates we are calculating approach 
the error in our numerical scheme, hence the relative er- 
ror becomes an issue. In fact after the 960th grid point 
the values we obtained for the derivatives of the Green's 
functions are badly behaved. So we ignore these data 
points and just simply interpolate between the 960th grid 
point and the values we will obtain for the derivatives on 
the horizon in Paper 1. 



C. Regularity of the Hartle Hawking State 

In Paper 1 we reduced the constraints on the stress 
tensor for the regularity of the state to requiring that 
the components of (Tj,^)ren possesses a Taylor series to 
the first order about the horizons of the black hole space- 
time of interest. As we have calculated both {Tj,^)ren and 
(Tj *)ren On the entire exterior region of a lukewarm black 
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hole, we may numerically calculate the radial derivative 
of each, either by using a finite difference method, or 
by fitting an interpolating function to the data points 
and differentiating. In Fig. [7] we plot the derivatives of 
{Tj.^)ren, (7t*)ren and (T'g^)reTi in the region of the event 
horizon. It is clear that both numerical derivatives are 
finite there. If this is the case then they must both pos- 
sess a Taylor series, at least to 0{r — rg), about both 
horizons. Therefore we may conclude, with correspond- 
ing confidence, that the equivalent of the Hartle-Hawking 
state for a lukewarm black hole configuration, is regular 
on both the event and cosmological horizons. 




FIG. 7: Radial derivatives of (T^ {Tg'')ren and {T,.'')re 

in the vicinity of the event horizon 




_ ('■- >■/.) 



0.7745 i 



CONCLUSIONS 



In this paper, we have, using the Hadamard renormal- 
ization procedure, developed a new prescription for the 
calculation of the {T^i,)ren for a spherically symmetric 
space-time possessing a thermal state. In this new ap- 
proach, each derivative of the Green's function is renor- 
malized individually, after which they may be combined 
to form the {T^^)ren- This method is equivalent to the 
method of Anderson, Hiscock and Samuel, but has the 
advantage of being more direct. Our method avoids the 
separate identification of divergences via the WKB ap- 
proximation and provides much greater control in debug- 
ging numerical calculations. 

The application of this new method naturally splits the 
{T^u)ren vnto two compouents; an analytic contribution 
which is formed from a combination of closed form ex- 
pressions, and a numerical component. We showed that 
the analytical component is a conserved quantity for a 
general spherically symmetric space-time. We then cal- 
culated, using Mathematica, the numerical component 
for a field in Hartle-Hawking state on the exterior re- 
gion (excluding the immediate vicinity of the horizons) 
of a lukewarm black hole, and demonstrated that it too 
is conserved. To facilitate this calculation we developed 
a method of performing the required summations of the 
WKB approximation of arbitrary order. 

Finally combining the results of this paper with the 
horizon calculations contained in Paper 1, we have plot- 
ted the components of (T^iy)ren on the exterior region of 
a lukewarm black hole. The combination of the results 
from both papers provides us with compelling numerical 
evidence that the equivalent of the Hartle-Hawking state 
for the lukewarm black hole is regular on both the event 
and cosmological horizons. 



FIG. 8: Radial derivatives of (T^ {Te'^)ren and {T/)re 

in the vicinity of the cosmological horizon 



Appendix A: Renormalization Subtraction Terms 

Here we list the renomalization subtraction terms for 
temporal separation. We recall here that if we are dealing 
with a space-time which has a constant Ricci scalar, such 
as the lukewarm case, then we can define an effective field 
mass, rh — + — l/6)i?. 
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g4j g2 I440r4 
-4r/ (2/' + r (rV" + 3r/"' - /")) + 4/2 - 4 (45mV4 + l) } In (-A/eV2) 

+ {3rV'' - 12rV/" (/" + Sm^) + 24r/3 (2/' + r (r^"" + 3r/"' - /")) 

+4/2 (r^ (SrV (/" - 30m2) - 4/'^ + r/' (9r/"' + 28/" - ISOrri^)) + 270TO'^r^ + 6) + 24/^} 



(A2) 



2/' 2/' 

OTT {g LjEsmg-tr' \ = g u~ i - / _ - -i — - 



^^^{g^* GEsing-rt'} — — 



,3/2- 

(/" + 60to2) + 4r2/'2 



(A3) 



e^P Sr^e^P 1440r4/ 

-2r3/' (r/'" + 2/" - GOm^) + 4r/ (r (r/'" + /") - 2/') + 4/^ - 4 (45m'*r'* + l) } In (-A/eV2) 

- ^^^^ {^Ir^r - ^r^ff" (26/" + Ihm') + 4/^ (lOrV" (/" + + r^/' (5r/"' - 24/" + BOm^) 

-90m V'' - 2) - 8r/^ (26/' + r {ir^f + 15r/"' - 25/")) + 8p} 



(A4) 



_pfn J (^pjii 2rf + 6m2r2 + 2) - 2p In {-Xfe^/2) 



6r^€^P 1440r4 
27- V (r/"' + 2/" + 60m2) + l20m^P - 4 - 2r/ (4/' + r {r^ f" + 2rf" - 2/" + 60m2)) + 4/^ + ISOmV^} 



„3 fi till 
1 

2880r4/2 



{llrV* - 8r2/3 (r (9/"' + r/"") + 16/") - 2r^ f {llP p + lOr/' + SOttiV^ + lO) + 1 
+4/2 (r (5r/" (r^/" + Gm^r^ + 2) - ISr/'^. + /' (5rV" + Wr^/" + 12077i2r2 + 20)) - 90m V"* - 2) } , (A5) 

+4r/ (2/' + r {P p" + 3r/"' - /")) - 4f + 180mV4 + 4} In (-A/e2/2) 

+ 2880r*p {-27^'/'" + 12rV/" (4/" + 25m2) + 4p (3 {Pp (/" - 30™^) + 24rf (2/' + r (r^/"" + 3r/"' - /")) 
-24/4 ^ 9Q^4^4 _^ 2) + 26r2/'2 - 2pp (3r/"' + /" + gOm^)) } , (A7) 



8Tr^{g^^GEsing;ee} = 



^2 J/2 _ J (^pjll _ 4^ J/ ^ g^jj2^2 + 2) + 2p 
gj,2g2 J2 



1 



1440r4 



{rV" + 4r2/'2 - 2pp {rp" + 2/" + 60m2) - 2rf {Ap + r (r^/"". + 2r/"' - 2/" + GOm^)) 



+4/2 + 180mV4 + UOfn^P - 4} In (-A/e2/2) + 2ggQj,^f2 {-H^V^ + 2r2//'2 (l7r2/" - 20r/' + 30m2r2 + lO) 

+4/2 {--dPp' {Pp' + 6m2r2 + 2) + 13r2/'2 + rf {-bP f" + 16r2/" + 60m2r2 - 20) + 90m V"* + 2) 

+Srp (30/ + r {Pp" - 6r/"' - 14/")) - 8/^} , (A8) 



{g'^'^GEsing-<j}<ji} — {g^^GEsingfie}- 
I 



(A9) 



Note we choose not to include the expression for 
g^^GEsing;rr as it wiU not be required. 



Appendix B: WKB Integrals 



In this appendix we derive the expressions (3.28) and 
([3^. 
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To derive Eq. (3.281 we firstly we consider the integral 
along Ca- We may isolate the divergent behaviour of the 
integral by initially subtracting off the Talyor series of 
gi about q = 1 from gi to the order which will render 



the integrand integrable at g = 1 and then adding back 
on the integral of this series over the denominator (1 — 
pollowing this procedure we obtain the following 
expression for the integral along Ca' 



n 



i-iy 
(1_ )«+l/2 



(Bl) 



For the integral along Cb we transform integration variable from g to in the following manner; g — )■ 1 + ee'^, then 
our integral becomes 



7^ 



fi/2 (_ge»e)'+i/2 



-d0 



(B2) 



We may then expand g,i{l + ee'^) about e = to obtain an expression which is now amenable to integration. This 
procedure gives the following form for the integral along C2: 



(-lp6^-+25^(l) 



a \ ^ 



which, together with Eqn (|B1 ), allows us to obtain 



r/2 ' - 



(1-9^+2 



(B3) 



(B4) 



As each term in Eq. ( |B4[ ) is manifestly finite as e — > we make take this limit leaving us with Eq. 3.28| 
For the derivation of Eq. (3.29) it is useful to introduce a new function 

j=o ■'' 

Then we have that G(l) — G"(l) = .. = G'^*^^''(l) = and so integration by parts yields 

1 G(q) _ G(0) G'(0) , G(*-i)(0) {-ly 



+ ■ ■■ + {-!)' 



(B5) 



1 GW(<7) 



(,_!)(,_!). (l-q)i/2 

(B6) 



Now using the definition of G{q) and Eq. ( 3.28[ ) we have 
^ a ( g,(0) ff,:(0) 



+ ... + (-ir 



+3.(1) 

+3:(i) 
+3r(i) 



* 2 * 2 



1 1 



1 1 



,_|2! *~i2! (^-i)(^-2)l! (, _ 1 )(, _ |)(, _ |) 



+(-l)'-'5r'(l) 



11 11 

I(T^ + z-i(z-l)! 
JO 



(z-i)(z-i)....(i)7o (l-g)^/^ 



(B7) 
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All the quantities in square brackets can be shown to vanish. 
Eq. (|3^ 



Therefore we arrive at the desired result, namely 
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